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^ . Abstract: Observations of pulsar timing provide strong constraints on scalar-tensor the- v — J 

QQ \ ories of gravity, but these constraints are traditionally quoted as limits on the microscopic ^J 

•/^ ■ parameters (like the Brans-Dicke coupling, for example) that govern the strength of scalar- ■v_^ 

matter couplings at the particle level in particular models. Here we present fits to timing ^ — >^ 

^— s ■ data for several pulsars directly in terms of the phenomenological couplings (masses, scalar 

^-H . charges, moment of inertia sensitivities and so on) of the stars involved, rather than to ^ — ^ 

. . . the more microscopic parameters of a specific model. For instance, for the double pulsar C_!) 

. ^ \ PSR J0737-3039A/B we find at the 68% confidence level that the masses are bounded by 

/\ \ 1.28 < mj^jmQ < 1.34 and 1.19 < rrig/mQ < 1.25, while the scalar-charge to mass ratios 

d I satisfy IoaI < 0.21, |as| < 0.21 and |as — OaI < 0.002. These constraints are independent 

of the details of the scalar tensor model involved, and of assumptions about the stellar 
equations of state. Our fits can be used to constrain a broad class of scalar tensor theories 
by computing the fit quantities as functions of the microscopic parameters in any particu- 
lar model. For the Brans-Dicke and quasi-Brans-Dicke models, the constraints obtained in 
this manner are consistent with those quoted in the literature. 
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1. Introduction 



o 
o 



o 
o 

Although General Relativity (GR) has many applications in astrophysics and cosmology, ^^ 

incomplete knowledge about the gravitating bodies involved prevents using most of these ^—-'^ 

as tests of the theory itself for the vast majority of systems outside of the solar system. ( — ^ 

Binary pulsars provide the rare exception to this rule, due to the great precision with which 

timing measurements allow the properties of their orbits to be inferred [Q, ^, ^. These 

orbital properties are parameterized by the inferred values of the Keplerian parameters 

(like semi-major axis, a, and eccentricity, e) that define the characteristics of the Newtonian 

orbit of the pulsar and its partner, together with its orientation in space. They are also 

parameterized by a suite of post-Keplerian (PK) parameters that describe observable slow, 

secular orbital changes over time, as well as relativistic time delays in the propagation of 

radio signals emitted by the pulsar. 

1.1 The classic constraints 

GR predicts the values for these post-Keplerian parameters in terms of the underlying 
Keplerian parameters and the masses m^.g of the two orbiting bodies. To test these 
predictions, bands are drawn in the tti^ — mg plane, of the form 

C*^m^,m,) = e°"^±A, (1.1) 

where ^^^ denotes the theoretically predicted value of a PK parameter, ^°^*^ denotes its 
observed value and A denotes the observational error. The validity of GR requires all such 
bands to have non-empty intersection, providing a substantive test provided at least three 



PK parameters can be measured. The area of mutual overlap then gives the GR-inferred 
masses of the stars, to within some tolerance. 

Thus far, general relativity has passed these stringent tests, using a number of bi- 
nary pulsars. For two of these — PSR 61534+12 |^ ^ and the double binary, PSR 
J0737-3039A/B |^, |^ |8| — the tests are particularly redundant since it is possible to infer 
observational values for five independent PK parameters. It is also possible to infer the 
spin-orbit precession frequency [^, ^, which can be considered as a sixth PK parameter. 
Moreover, for the double binary, PSR J0737-3039A/B, both stars are pulsars, and a mea- 
surement of the ratio of the semi-major axes yields another constraint on the masses. The 
masses are constrained quite accurately 



?7i^ = 1.3381(7) m0 and 



mr 



1.2489(7) m© (GR) 



(1.2) 



The success of GR in describing pulsar orbits also constrains alternative theories of 
gravity; requiring their predictions to agree with GR to within present errors. Prominent 
among these alternatives are scalar-tensor theories H, [10], for which long-range gravita- 
tional forces are mediated by both the metric, g^^, and a very light scalar, (/), described by 

the action 
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{R^u + df,(f) dy 
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(1.3) 



Here k^ = ^ttG/c^ = h/{M'^c^) denotes the gravitational coupling, while Sm = Sm[g^j_y, (j), ^j] 
denotes the matter action, and controls how (j) and g^^u couple to the various other 'matter' 
fields, ^j, which we take in what follows to have the form 



o 
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o 
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Sm = Sm[A^i(l))gfMU,'^i 



(1.4) 



This form of coupling has several motivations. First, it is favored by strong observational 
constraints [O, 12 1 on violations of the weak equivalence principle, which are evaded by 
actions of the form of eq. (|l.4|). Second, it is also the kind of theory that actually describes 
the low-energy limit of certain types of extra-dimensional theories ||l3|. 

The predictions of theories of this type have have been compared in detail with binary 



pulsar data [14, |15|, with the coupling function assumed to have the particular quasi-Brans 
Dicke form, ^ 



A^cp) = exp 



«s </* + Y ' 



for which the effective coupling of scalars to matter turns out to have the strength 



«(</>) : = 



dlnA 



as + bs(j). 



(1.5) 



(1.6) 
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This form is motivated by the idea that the field (j) does not vary appreciably in any 
particular system of interest, and so a{(j)) is approximately constant. 



^Brans Dicke theory corresponds to the specific choice bs — 0. Our notation follows that of ||l 
differs in minor ways from that used elsewhere in the literature ||lOt|, as described in detail in table 
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Our Notation 


Notation of 


m 


Meaning 


as 


a 




Microscopic scalar-matter coupling constant 


bs 


P 




Microscopic scalar-matter coupling constant 


a{cp) 


«(</>) 




Microscopic scalar-matter coupling function 


aA,B 


(y-A,B 




Effective scalar-matter coupling of a star 


bA,B 


Pa,b 




Effective scalar-matter coupling of a star 


Ka,b 


I^A,B 




Sensitivity of a star's moment of inertia to scalar field 


Qa,b 


^A,B 




Scalar charge of a star 



Table 1: Table of Notation. 

Comparison with solar system and pulsar data is found to constrain the microscopic 
couplings for this theory, a^ and bg, to be consistent with zero, with a^ strongly constrained 
from solar-system data and hg restricted by pulsars. The constraints on these microscopic 
couplings are usually presented as an exclusion plot in the bs-dg plane, and in particular it 
was found that [0, 15 1: 
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-4.5. 



+ bs 






< 3.4 • 10" 



(1.7) 



where (f)^ is the value of the scalar field asymptotically far away from the solar system, 
which is conventionally taken to be zero. 

1.2 A more model-independent approach 

There are two related drawbacks to traditional comparisons between scalar-tensor theory 
and observations. First, because limits like eq. (1.7) are quoted directly for the microscopic 
couplings, as and 6s, a completely new analysis is required for each new assumed functional 
form for the coupling function A{(j)). Second, as discussed below (and remarked on by the 
original authors), these bounds are subject to uncertainties that are hard to quantify, due 
to limits of our understanding of the nuclear equation of state that applies within the 
pulsars. 

One approach towards robustness that has been taken in the literature is to generalize 
the form taken for the phenomenological lagrangian describing the two-body interactions 
of the gravitating objects, restricting attention to theories of gravity whose predictions for 
the orbital dynamics may be derived from a boost-invariant Lagrangian (at least to the first 
post-Newtonian (IPN) order). Will and Damour & Taylor have shown that such theories 
may be characterized by a set of body-dependent phenomenological parameters |12, @], of 
which there are five — m^, tti^, Q, e, S, — in the most general Lagrangian |12, g]: 



,.z!+« + i,a_3.r+«((.+„)^=+.(N.v)^ 



R 



8c2 



2Rc^ 
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where R and V are the relative position and velocity vectors, N = R/i?, and Q, uia, tub, 
e, and ^ are phenomenological parameters, and M = niA + ruB, v = mArriB/M'^. Wex 
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and Kramer have used the observed values of the PK parameters for the double pulsar to 
constrain these WDT (Will-Damour- Taylor) parameters ^. 

However, a drawback of this approach is that it does not describe radiation effects, 
and in particular the PK parameter Pb (which describes the shortening of the orbital pe- 
riod due to emission of gravitational radiation). This is because P^ cannot be expressed 
in terms of the WDT parameters. This limitation exists because radiative effects come in 
at higher order in powers of V/c (at 1.5PN order for dipole emission and 2.5PN order for 
quadrupole emission), whereas the WDT parameters only describe IPN effects. Extend- 
ing the phenomenological parametrization to higher PN orders introduces too many new 
parameters for the data to usefully constrain, however. Thus, in order to use all PK pa- 
rameters, including P^, to obtain interesting constraints on alternative theories of gravity, 
it is necessary to further restrict the class of theories that one considers. 

In the present paper we take a complementary approach to confronting pulsar observa- 
tions with scalar-tensor models, based on the observation that the scalar-tensor predictions 
for the PK parameters depend only on a relatively small number of macroscopic quantities C ^ 

that characterize how the two stars couple to the scalar field. There turn out to be seven ( — ^ 






of these, (defined in detail in section 2^): mA,B, cla,b, bA,B and /ca- These generalize the 

two masses — tjia and itib — that suffice to make predictions within General Relativity, to 

include four new quantities — 0^,5 and 6a,s — that characterize the strength with which '■^^ 

the scalar field couples to the pulsar and its orbital partner, plus one variable, kA, related C_-) 

to the pulsar's moment of inertia. C ^ 

The key point is that model-dependent complications, like the detailed form of A{(f)) ^1 

and knowledge of the nuclear equation of state, enter only into the predictions for the ^_^ 

quantities niA^B, cla,b, ^a.b and Ua as functions of the microscopic couplings in a particular 
scalar-tensor theory. But these complications do not enter at all into the formulae that 
express how m,A,s through A;^ determine the observed PK parameters. v — J 

This makes it useful to phrase the confrontation between theory and experiment in two C_) 

steps: first use the observations to constrain the quantities niA^B through ^4 once and for 
all in a model-independent way; then compute these quantities within specific scalar-tensor 
models as functions of the underlying model parameters that define the function A[(j)). 

It is the goal of this paper to perform the first — and model-independent — one of 
these steps. At first sight this might seem to be impossible to do, since it appears to 
involve constraining more quantities than the five observable PK parameters^. However, 
it turns out that the present data nonetheless allow useful constraints to be achieved, for 
two reasons. First, constraints are possible because the dependence of the PK parameters 
on two of the parameters, 6a,b, is very weak; thus effectively reducing the number of free 
variables from seven down to five. 

Second, one combination of Keplerian parameters — the ratio of projected semi-major 
axes of the two orbits, R = Xb/xa (see section |2?T| for details of notation) — can sometimes 
also be measured. For instance for the double pulsar observations give R = 1.0714(11). This 



^ Since the theoretical prediction for the spin-orbit precession frequency has not yet been calculated in 
scalar-tensor gravity, we do not include it in the present analysis. 



measurement is useful because the theoretical prediction for this quantity in all Lorentz- 
invariant theories of gravity is Rth = mA/mg + 0(l/c^) and so is completely independent 
of the scalar couplings aA,B and /ca |jl5| (see also section p.l[ ). When this ratio is measurable 
the number of observational constraints to be satisfied rises from five to six. 

In the end we find that interesting model-independent constraints on m^ ^ and a^^B 
are possible. The great virtue of these bounds is that they directly constrain the stellar 
parameters on which the PK parameters depend, independent of any assumptions about 
the function A{(p) and the nuclear equation of state. 

Once the observational constraints on these quantities are known, they can be used 
to constrain the microscopic parameters for any choice of A{(f)) or equation of state as a (~ , 

separate step. When we do so for the quadratic model A{(j)) = exp(as0 + 6s</>^/2), we find 
agreement with earlier work, which already rules out a large part of the most interesting 



The resulting bounds are then subject to uncertainties in the equation of state, which are 
hard to quantify. 

In the end, the phenomenological analysis to which we are led in this paper is similar 
in spirit to that carried out by Wex and Kramer ^. The difference is that we specialize 
to a more restricted class of theories of gravity, for which the energy loss due to emission 
of gravitational radiation has been calculated and takes a relatively simple form. Conse- 
quently, unlike Wex and Kramer, we are able to make use of all observed PK parameters, 
including Pb, when obtaining constraints. 

We thus urge observers to express their results in this more model-independent way, 
which potentially can then be used by theorists to constrain a great variety of specific 
models. Its independence of the internal structure of the objects involved also shows that 
some of their physical properties - like the masses mA,B ~ can be inferred quite robustly, 
without making assumptions about which theory of gravity actually applies in Nature. 






region (that of spontaneous scalarization jTj, 15, p^]). 

It might come as a surprise to the reader that this second step might depend on the 
details of the neutron equation of state, since it is an important feature of general relativity C ^ 

that predictions for the PK parameters depend only on the two masses, itia^b, of the stars, ^~^ 

and on none of their other properties. This happy property is called the 'principle of 
effacement' of internal structure [T^ , and it is this feature that allows a precise prediction 
of all PK parameters in GR using only the masses, without need for detailed knowledge 
about the star's structure. C — J 

In scalar-tensor gravity (as in most other alternative theories of gravity) the principle C ^ 

of effacement does not hold. The prediction for the PK parameters actually depends in ^1 

principle on all seven of the quantities which characterize the internal gravitational fields: ^_^ 

iTT'A^B, (^A.B, bA,B, ^A- For uow, it is important to note that these are all a-priori independent. 
Once a particular equation of state is specified, then the equations of stellar structure can 
be solved, and on a given branch^ of stellar configurations, the quantities a, b and k for ^« — J 

each star can be expressed in terms of its mass m and the underlying parameters — like C !) 

Og and hg — that define the scalar-tensor model. This is the approach taken in [14, 15|. 



■^In scalar-tensor theories there may be multiple branches of stellar configurations Ifj 
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2. Formalism 

In this section we collect expressions for the Keplerian and Post-Keplerian parameters in 
scalar-tensor theories, following the results of [§, y] and references therein. We do so both 
to establish notation and to provide context for the bounds obtained in the next section. 
Experts and readers in a hurry should feel free to skip this part completely. 

2.1 Keplerian and Post-Keplerian Parameters 

First, a reminder of how Keplerian orbits are described, followed by the post-Keplerian 
parameters that describe slow secular changes to the Keplerian parameters. _ 

M 

Orbital description ^~\~* 

The non-relativistic gravitational two-body problem famously predicts bound orbits to PJ 

be ellipses. More specifically, in the center of mass frame the relative position vector, r = hT~J 

r4 — "Tg, sweeps out the trajectory of an ellipse whose shape is described by two parameters: 

the semi-major axis o, and the eccentricity e. The positions r^ and fg also separately trace 

out ellipses, with semi-major axes satisfying a^/a = mg/M and Og/a = mjx/M, where 

M = tua + rng is the system's total mass. -^"^ 

The time taken to traverse this orbit is given by the orbital period Pf,, and is re- [\^ 

lated to the semi-major axis by Newton's modification of Kepler's Third Law:^ Ph = < — ^ 

27ra^'^(^M)~^'^. For timing measurements a reference time, Tq, is also needed to specify 
the time of passage through periastron. 

The orientation of the orbital ellipse with respect to a reference triad is specified by 
three angles. In celestial mechanics these angles are conventionally taken to be the longitude 
of ascending node, r/; the orbital inclination, i, of the orbit relative to the plane of the sky; C ^ 

and the argument of periastron, to. r — n 

Standard techniques allow the positions and momenta of the two bodies to be inferred 
from the variables (a, e, r\, i, to, Tq) by means of the action-angle formalism [ p!s[| . 

Deviations from the Newtonian two-body problem cause the shape and orientation of 
the Keplerian orbits to change with time. In practice this change is slow enough to be 
regarded as a small secular evolution in each of the Keplerian parameters. Two of these 
have been accurately measured for several binary pulsars: the decrease in the orbital period 
due to emission of gravitational radiation, Pi^; and the precession of the orbital periastron, 
u. 

Pulsar timing 

Pulsars emit electromagnetic signals at regular intervals that are measured on Earth by 
radio telescopes. The problem of relating the time of emission T^ at the pulsar to the time 
of arrival Tq on Earth is usually split into two parts. The first part involves relating Tg to 
the time of arrival ta at the solar-system barycenter, neglecting the time delay due to the 
solar gravitational field. The second part involves relating ta to r^, which depends only 



o 



o 



*5 is the effective gravitational coupling constant between the two bodies. In scalar-tensor gravity, 



on the motion of Earth with respect to the solar system barycenter (and is not considered 
here) . 

The solution to the first problem is the pulsar timing formula, and is conventionally 
written 

Z?ta = Te + A^ + As + A5 + 0(1/C^) . (2.1) 

These terms have the following origin. 

• The Doppler factor D describes the time dilation caused by the motion of the solar 
system relative to the binary pulsar. 

• The term A^ is called the Romer time delay. It is ©(l/c), and is due to the depen- ( | 

dence of the light path on the position of the pulsar. The variable that sets the time \~V^, 

scale of the Romer delay is the projection of the semi-major axis along the line of \~r\ 
sight, measured in units of time. It is conventionally denoted Xa.b = i^.s sini/c, and 
is called the light crossing time. 



2.2 PK Parameters in Scalar- Tensor Gravity 

When solving the problem of stellar structure in scalar-tensor gravity, it is necessary to 
specify the boundary condition for the scalar field asymptotically far away from the star, 
0OO. All of the properties of the star, such as its mass, therefore depend implicitly on (ft^o- 
(Because these properties can be multiple- valued functions of (/)oo, for some purposes it 
can be useful instead to follow the dependence on the value of the scalar field at the star's 
centre [|16| .) 

Although the asymptotic field is simple to specify for isolated stars, it is a more compli- 
cated concept for an orbiting binary system. For isolated stars, the scalar field sourced by 
the star becomes small far from its position, leading to the generic weak-field large-distance 
form 

^..(t,0 = ^. + ^3|^ + ..., (2.2) 



o 
o 



The term A^ is called the Einstein time delay. It is ©(l/c^), and is due to the 
time dilation caused by the motion of the pulsar A and the gravitational field of the 

companion B. The variable that sets the time scale of the Einstein delay is called 7, """^ 

and is given by equation ( 2.1C| ), below. LxJ 

The term Ag is called the Shapiro time delay. It is 0{l/(?)^ and is caused by the 

effect of curved space on the propagation of light. It can be thought of as the first ^ — ^ 

relativistic correction to A^. The variable that sets the time scale of the Shapiro delay ^^ 

is called r (the range of the Shapiro delay), and is given by equation ( ^.11 ) below. ■**— ^ 



It is proportional to Goos, the gravitational coupling between the companion B and /- — n 

a photon, whereas all the other PK parameters depend on Gab, the gravitational 
coupling between the two orbiting bodies. It is intuitively clear that the importance 
of the Shapiro effect depends strongly on the inclination angle of the orbit. The 
variable which parametrizes this dependence is called s (the shape of the Shapiro 
delay), and is given by equation ( 2.12| ). 



o 
o 



where Qa,b defines the scalar 'charges' of the two stars, and rA,B{t) are their trajectories.^ 
We here make the choice that the scalar field approaches $00 asymptotically far away from 
both stars of the binary system. It is convenient to choose units of length so that <^oo = 0. 
Now consider a binary system of stars A and B, and assume that the separation 
between the two stars is large enough to justify the near-Newtonian weak-field limit. In 
this regime the fields sourced by the two stars can be superposed so the total scalar field is 
$ = $^ -|- $B- Consider now the region much closer to one of the two stars, and let 4>a and 
4>B denote the scalar fields of each star in this 'internal' regime. The boundary conditions 
for this 'internal' field therefore are (approximately) 

These expressions show that (i;^a,b)oo are small whenever Gm/rcP' ~ {v/c)'^ <C 1, where r I J. I 

is the separation between the two stars, and v is their orbital velocity. h— J 

Now, the couplings that are relevant for computing post-Keplerian quantities govern 
how the mass of star A depends on the boundary condition (<^a)oo- For instance, expanding 
777,4 in a power series about (</>a)oo = 0, defines the coefficients a a and 6a: 

"T'a[('/'a)oo] = m^ 



OooW^$s(i,rA(t)), {4>B)oo{t)^^A{t,rs{t)). (2.3) 



A;a [(</*/ 



d{(t>A 
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l + aA(</)A)oc + ^(6A + (aA)')(0A)L + --- , (2.4) l>0 

o 

and an expansion of tub about {4)b)oo = similarly defines Og and 6^. We may define the (^~-^ 

scalar coupling functions, Oa.b and hA,B-, by 1 

\(A. \ 1 - d\0gmA[{(t)A)oo\ ,UfA\^ for.^ ' 

fflA ((PAJOO = ^7—^^ = aA+OA[(PA)oo + ■■■ , (2.5) 

0{(Pa)oo (^ 

6a[(0a)oo] = r.,: .2 = 6a + ..., (2.6) L^ 

(^WaIoo ^,— >. 

with similar definitions for aB[{(j)B)oo] and 6s[(i?^s)oo]- It is a general property of scalar- 
tensor systems that couplings defined in this way agree with those defined from A((f)) using 
eq. ( |1.6D in the limit of weakly coupled non-relativistic systems. 

Pulsars rotate, and the frequency of pulsation is given by the rotational frequency 
$7 = J/I, where J is the 'spin' angular momentum of the pulsar, and I is its moment of 
inertia. I can be found by solving the equations of stellar structure, and just like the mass, 
it depends on the boundary conditions for the scalar field: 

/a[(0a)oo] = /a [1 - kAiMoo + ...], (2.7) 

leading us to define 

_ d log I A[{(t> a) 00] , , ,„ „x 



''Knowledge of the stellar equation of state allows Qa.b to be computed as functions of the corresponding 
mass, VTiA.B, and the asymptotic value of the scalar field. 



The predictions for the PK parameters in Scalar- Tensor gravity turn out to be deter- 
mined by the Keplerian parameters, together with the coefficients m,^ a, a^^g, b^g and fc^ 
ra. Exphcitly, 
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GabMtA 
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1 + a^Os 
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(Xb(1 + OAais) + 1 + /cAflf 



(2.9) 

(2.10) 
(2.11) 

(2.12) 



where G^b = G{1 + a^aB) is the total (graviton plus scalar) weak-field coupling between 
the two bodies, and G^b = ^(1 + a^^aB) is the coupling between the companion B and 
a non-compact body in the vicinity of the binary pulsar. The quantity a^^ is given by 
the value of the scalar-matter coupling a{(j)) (defined in equation ( |1.6| )) asymptotically far 
away from the binary pulsar. As before M = uia + m-s is the total mass, while n = 27r/Pfe 
is the orbital angular frequency and X^^b = "Ia,s/-^- 

The expression for the decay of the orbital period is similarly given by 



Pb - ^b 
where the different contributions are 

SttXaXb 



pmon I pdip Aquad pkin , pgal 



(2.13) 
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,Mny/^ (1 + 73eV24 + 37e^/96) 
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(2.14) 
(2.15) 

(2.16) 



where 'mon', 'dip', and 'quad' denote monopole, dipole, and quadrupole radiation, respec- 
tively. In the limit of general relativity — i.e. when aA,B — ^ and 6.4. s — )• — the 
monopole and dipole contributions vanish. Note that the monopole term is of order 1/c^, 
because the total scalar charge of the binary system is constant in time. 

The last two terms in equation ( 2.13| ) arise due to the relative motion between the 
binary pulsar and the solar system [20|. The kinetic contribution to Pb is given by 






kin 



Vj,/cd, 



(2.17) 



where v is the velocity of the binary pulsar relative to the solar system, and T denotes the 
component transverse to the line of sight, and d is the distance between the binary pulsar 
and the solar system. 
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The galactic contribution to Pi, is given by 

Pr' = ««/c, (2.18) 

where a is the acceleration of the binary pulsar relative to the solar system, and R denotes 
the component along the line of sight. We have g, ^ 



1 I gal 

1 + (aio"^)2 



gpsr-gss, (2.19) 



where ^psr is the acceleration of the binary pulsar relative to the galactic centre, as pre- 
dicted by a Newtonian galactic model, and ^ss is the corresponding quantity for the solar C i 
system. The quantity af^ is the value of the scalar-matter coupling function a{(p) (defined I— i-i 
in equation (|1.6D) asymptotically far away from the galaxy, and 



Opsr = XaUa + Xgag (2.20) h^ 

is the charge-to-mass ratio of the binary pulsar system as a whole. ^« — -' 

A simple galactic model may be used to relate a^ to a^'^ , and so in principle the PK C !) 

parameters depend on a^ (and the galactic model) in addition to the seven quantities ^_,^^ 
iTT-A.B, cla.bi ^a,b and k^. However, it turns out in practice that the contributions of a^ 
are much too small to be measurable, so we henceforth set a^ = a^ = 0. This ensures 
that GR is recovered asymptotically far away from the binary pulsar, and asymptotically 

far away from the galaxy. C J 

Note that equations ( p.9[ )-( 2.13 ) are invariant under ^J 



O 



flA -^ -a-A , clb -^ -as , a^' -^ -a^' , af^' -^ -a|^^ and k^ -^ -kA , (2.21) 

o 

which corresponds to switching the sign of the scalar field. 

3. Constraints ^>— ^ 

3.1 Statistics 

In this section, we use the method of least squares [£^] to compare the observed values of the 
PK parameters to the predictions of scalar-tensor gravity, thereby obtaining constraints on 
the phenomenological stellar parameters {m-A.s, aA.B,^^,^, ^a}- For brevity, we will denote 
these phenomenological stellar parameters by {PjIJ^^. 

For a given pulsar of interest, let {Ci}iLi run over as many of the quantities i\, cj, 7, r, 
s and R as are measured. Assume that the measurement process for ^j can be described by 
a normal (Gaussian) distribution with standard deviation Aj, and assume that correlations 
between the different ^i can be neglected. 

If we know a priori that the pulsar timing is correctly described by the phenomenolog- 
ical stellar parameters {Fj} within scalar-tensor gravity, then the probability of measuring 
any given set of values of {^j} is given by 

P(ft,...,&|r....,r,) = ^^-^,,l^^, (3,1) 
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where 

" "Ci-e(ri,...,r7 



^' = E 



N r> .th/ 



tth 



(3.2) 



where Q are the theoretically-predicted values in scalar-tensor gravity. By means of Bayes' 
theorem, this probabilistic statement may be turned around: if we know a priori that a 
given set of values of {^j} have been observed, then the probability that the pulsar is 
described by the phenomenological parameters {rj} in scalar-tensor gravity is given by 

p(ri,...,r7|ei,...,eiv)= pi^'''"'J'i p(6,---,g7v|ri,...,r7). (3.3) 

Assume that we have no prior information about {^j} or {rj}. Then, combining equations i— I— i 

dp) and (|3t) yields \jj 

P(T^,...,Tr\^,,■■■,^N)r-e-^'/\ (3.4) hj^) 

In principle, equation ( |3.4D may be calculated numerically, and integrated over the seven- C_) 

dimensional parameter space to find the constraints of interest. For instance, the mean (^ ^ 

value and variance of Fj are given by 



O 

, /d^r(r.-r.)^e-xV2 ^ 



and 



respectively. In practice, it is very difficult to calculate ( |3.5| ) and ( |3.6| ) directly. Therefore, C ^ 

we resort to approximation methods. Assume that x^ has a global minimum at {T*}, and /- — ^ 

approximate it near this minimum by a quadratic form: 



X'(ri, . . . , Ty) = Xmin + E (^^ " rn^.^CT, - r,*) . (3.7) 



Substituting equation (3.7) into equations (p.4p, (p?5|), and ( pTq ) shows that the Pj are 
normally distributed, with means 

f. = ^^ (3.8) 

and variances 

^r. = 7^ • (3.9) 

The off-diagonal components of C are related to the correlation coefficients between the 
different {TJ. 

3.2 Implementation 

We have implemented the statistical analysis described in the previous section for two 
pulsars, for which at least five PK parameters have been measured. For the double pulsar 
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o 



PSR J0737-3039, all of the quantities Pb, uj, 7, r, s, R have been measured, and so A^ = 6. 
For the binary pulsar PSR B1534+12, by contrast, R is not measured and so A^ = 5. 

The double pulsar PSRJ0737-3039A/B g 0, |] consists of two pulsars, A and B, 
that are bound in a relativistic orbit described by the Keplerian parameters summarized 
in table ^, and the Post-Keplerian parameters summarized in table ^ Pulsar A has a 
period of 23ms, and pulsar B has a much slower period of 2.8s. The Post-Keplerian timing 
parameters 7, r, s all pertain to pulsar A. The light crossing time x has been measured 
for both pulsars. For the binary pulsar PSR B1534-I-12 Q, the Keplerian parameters are 
summarized in table §, and the Post-Keplerian parameters are summarized in table ^. 



Symbol 


Meaning 


Value 


X \ 


Light crossing time of A 


1.415032(l)s 


Xg 


Light crossing time of B 


1.5161(16)s 


e 


Eccentricity 


0.0877775(9) 


Pb 


Orbital Period 


0.10225156248(5)days 


UJ 


Argument of periastron of A 


73.805(3)° 


To 


Time at periastron of A 


52870.0120589(6)MJD 



o 
o 



Table 2: Keplerian parameters for the Double Pulsar J0737-3039. 

For these two pulsars, we have calculated (^]^) numerically on a grid in the seven- 
dimensional parameter space, and found that x^ has a global minimum near Ua = clb = 0. 
This is not surprising, because these pulsars are very well described by general relativity. 
We have also found that the minimum value of \^ is very close to zero ~ Xmin ~ 10~^ for 



the double pulsar, and XZ. 



10 ^ for 1534-1-12. This is also not surprising, because we 



have more parameters than data points. 

Since it is hard to visualize the seven-dimensional parameter space, we will define the 
following functions in order to present the results of our numerical calculations: 



O 
O 

O 
O 
O 



Xm(m^,?ns) = min x ("T'A,mB,aA, Os, 6^, &s, ^a) , 

xliaA,aB) = min x^imA,mB,aA,aB,bA,bB,kA) ■ 

mA,mB,bA,bB,kA 



(3.10) 
(3.11) 



If we use the quadratic form approximation ( |3.7| ), then Xm ^-nd Xa ^^^ ^^^o quadratic 
forms, whose coefficients can be explicitly calculated in terms of the components of C p2|. 



Moreover, it can be shown that [22| 



minXa(«A + S,aB) = Xmin + "^V^",^ 



as 



min xliaA, a% + 6) = xLin + ^'^/^Ig 



aA 



Xmin 
Xmin 



aA 



(3.12) 
(3.13) 



and similar relations hold for Xm- Thus, the standard deviations may be found from contour 
plots of xl and Xm- 

For the double pulsar PSR J0737-3039A/B, contours of Xa are plotted in the left panel 
of figure |l]. This plot shows that a a is very close to a^ , with a precision that is more easily 
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Symbol 


Meaning 


Value 


Pb 


Time derivative of orbital period 


-1.252(17)10^^2 


Co 


Precession frequency of periastron 


16.89947(68)°yr-i 


7 


Relativistic timing parameter 


0.3856(26) 


r 


Range of Shapiro delay 


6.21(33);us 


s 


Shape of Shapiro delay 


0.99974(-39,+16) 



Table 3: Post-Keplerian Parameters for the Double Pulsar J0737-3039. The galactic and kinetic 
contributions to Pf, are negligible. 

seen in the right panel of the same figure. The robustness of this constraint can be simply 
understood from equation ( 2.15 ) - if |as — a^l gets too large, then this equation predicts too 
much dipole radiation, inconsistent with observations. Changing variables from (0^,0^) to 
(ca, as — Oa) in equations ( |3.12| )-( |3T^ ), we infer that 



0", 



OA 



0.21, 



0", 



aB-a-A 



0.002 . 



(3.14) 



while the mean values of both Ca and a^ — Oa are zero. Also, note that the plots in figure 
|l] are both symmetric under aA,B -^ —O'A.b, as expected from the symmetry (2.21). 



This same analysis may be carried out for the pulsar PSR B1534+12. The contour 
plots of Xa ^^^ shown in figure |2[ We find that 



0", 



a.A 



0.44, 



0", 



a-B-a-A 



0.004. 



(3.15) 



For the double pulsar PSR J0737-3039A/B, contours of Xm ^^^ plotted in the left panel 
of figure 1^. The masses lie very close to the line ttt-a = Rtub , in good agreement with the 
general prediction. The deviation from this line is shown in more detail in the right panel 
of the same figure, by plotting XbIUb — XaITIa versus ttia. The minimum value of Xm is 
close to the mass values inferred in GR — itia = 1.3381(7)r?T,0 and tub = 1.2489(7)m0 — 
at the top-right of the line in the left panel of figure ^, and in the center right of the right 
panel of the same figure. 

Note that the behaviour of Xm near its minimum value is very asymmetric. This 
asymmetry is caused by the contribution of equation ( p.l2[ ) to x^- At the minimum value 
of Xm^ ^^ have s ~ 1. The value of Xm increases very rapidly as we enter the region s > 1, 



o 
o 

o 
o 
^ 

o 
o 
o 



Symbol 


Meaning 


Value 


Xa 


Light crossing time of A 


3.729464(2)s 


e 


Eccentricity 


0.2736775(3) 


Pb 


Orbital Period 


0.420737299122(10)days 


UJ 


Argument of periastron of A 


274.57679(5)° 


To 


Time at periastron of A 


50260.92493075(4)MJD 



Table 4: Keplerian Parameters for PSR B1534+12. 
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Symbol 


Meaning 


Value 


n 


Time derivative of orbital period 


-0.137(3)10-12 


ib 


Precession frequency of periastron 


1.755789(9)°yr-i 


7 


Relativistic timing parameter 


2.070(2) 


r 


Range of Shapiro delay 


6.7(1.0);us 


s 


Shape of Shapiro delay 


0.975(7) 



Table 5: Post-Keplerian Parameters for PSR B1534+12. In the value reported for P;,, the galactic 
and kinetic contributions have been subtracted. 

whereas Xm increases much more slowly as we enter the region s < 1. Theoretically, s is 
the sine of the orbital inclination angle, i.e. s = sini and so s < 1. Strictly speaking, this 
constraint should have been imposed by a prior probability in equation ( |3.3| ). However, we 
see that in practice, for the double pulsar, this constraint is automatically enforced by the 
rapid growth of Xm- 

This strong asymmetry also shows that the quadratic form approximation to x^i sq. 
(3.7), is not a good one, and that strictly speaking, the use of equations ( 3.12| )- (|3.13| ) is 
not justified. However, we will still use the Xm — Xmin + ^ contour to estimate the allowed 
range for the masses at the 68% confidence level. We find 



(3.16) 



1.28 < niA/niQ < 1.34 , 1.19 < mg/mQ < 1.25 . 



Now for the pulsar PSR B1534+12, contours of Xm ^^^ plotted in figure 0. The minimum 
value is at niA ~ 1.23?tt,0, and mg ~ 1.29m0. The masses are not as highly correlated as in 
the double pulsar case, because the ratio R has not been measured. Also, the asymmetry 
of Xm isn't as pronounced as in the double pulsar case, because the PK parameter s has 
been measured much less accurately than for the double pulsar. Looking at the contour 
Xm = Xmin + 1) we find that at the 68% confidence level. 



o 
o 

o 
o 
^ 

o 
o 
o 



0.97 < mA/niQ < 1.28, 



1.15 < niB/mQ < 1.31 . 



(3.17) 



Note that in GR, the masses inferred from the combined measurements of w, 7, r, and s 
are ?tt,^ = 1.3332(lO)m0 and TUg = 1.3452(10)mQ [^, which correspond to the top-right 
corner of figure ^ and lie outside of the range ( |3.17| ). However, including Pf, significantly 
increases the errors of the GR-inferred masses. It is often difficult to obtain an accurate 
measurement of Pb, because this PK parameter receives galactic contributions. 
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Figure 1: Contour plot of x^ for the Double Pulsar J0737-3039 as a function of Qa and a^ 
(left panel), and a function of a^ — a a and a a (right panel). The value of x^ is minimized over 
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Figure 2: Contour plot of x^ for PSR B1534+12 in the a^ — a^ plane (left panel), and with Uj 
shown vs a A (right panel). The value of x^ is minimized over tti^, m^, bAjbg^kA- 



o 
o 

o 
o 
^ 

o 
o 
o 
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Figure 3: Contour plot of y^ for the Double Pulsar J0737-3039 in the m^ — vhb plane (left panel). 
The right panel plots the same information using a variable that emphasizes the accuracy of the 
test of the prediction R = mA/ras- The value of x^ is minimized over aA,aB,bA,bB,kA- Units on 
both axes are in solar masses. 



o 
o 

o 
o 
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Figure 4: Contour plot of x^ for PSR B1534+12. The value of x^ is minimized over 



aA,ambA,bB,kA- Units on both axes are in solar masses. 
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If we consider the quasi Brans-Dicke model with A{(j)) = exp(as(/' + 6s0^/2), and use 
the same relativistic polytrope models for neutron stars as in |l^] and |l^, then we find 
that the double-pulsar constraints ( p.14 ) and ( |3.16| ) together imply that bs > —5.5 for 



EOS II, and 5, > -4.6 for EOS A. For PSR B1534+12, the constraints (|3l5| ) and ( |3ll ) 



together imply that bg > —6.2 for EOS II, and bg > —5.4 for EOS A. These bounds on b. 



s 






are similar to those presented graphically in |15|. 

4. Conclusions 

In summary, we have shown that the existing data for two pulsars is constraining enough 

to place model-independent bounds directly on the stellar parameters that control the size 

of post-Keplerian eff'ects in scalar-tensor models. The virtue of these bounds is that they 

do not depend on the particular form for the function A{(j)) that defines which model is L J 

of interest, or on the details of the stellar equations of state. In particular we find that i"^ 

existing data impose strong and model-independent constraints on the relative size of the /- — >\ 

two scalar charges and masses. ^ — ^ 

These bounds can also be used to constrain particular models by computing in these 

models the masses and couplings as functions of the microscopic parameters. It is at this '""^ 

point that dependence on things like the stellar equation of state enters. L\-) 

We applied these methods in particular to the double pulsar J0737-3039A/B, and found C ^ 

that 1.28 < mA/rriQ < 1.34, 1.19 < mB/rriQ < 1.25, \aA,B\ < 0.21, and |as - a^| < 0.002, ^— ^ 

with 68% confidence, for all choices of scalar-matter coupling function, and for all nuclear - 
equations of state. 

A similar analysis of the pulsar PSR B1534-M2 yields 0.97 < m^/mo < 1.28, 1.15 < 

iriB/mQ < 1.31, Iua.bI < 0.44, and lag - a^l < 0.004. CD 

o 
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